Phase space topology of a switching current detector 
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I We examine in theory and by numerical simulation, the dynamic process of switching from a zero 

. voltage to a finite voltage state in a Josephson junction circuit. The theoretical model describes 

' small capacitance Josephson junctions which are overdamped at high frequencies, and can be applied 
to detection of the quantum state of a qubit circuit. We show that the speed and fidelity of the 

^ : readout are strongly influenced by the topology of the phase space attractors. The readout will be 

' close to optimal when choosing the circuit parameters so as to avoid having an unstable limiting 

pH ' cycle which separates the two basins of attraction. 

\ 

1 • One of the unique and exciting characteristics of quantum bits based on Josephson junction circuits (as opposed to 
^ ', quantum systems like atoms or nuclear spins) is the ability to engineer the system for optimal performance with regard 
^ to decoherence and readout efficiency. With Josephson junction circuits we have electrical access to the Hamiltonian 
c/2 ' of a quantum system described by two conjugate circuit variables. External current sources or voltage sources can 

be used to change potential and kinetic energies in the Hamiltonian of the quantum circuiti. This fact allows one to 
design a fast and reliable readout. 

The quantum dynamics of Josephson junction circuits is intensively studied by several groups around the world 
today. Many different circuits and measurement schemes have been proposed and implemented in experiments which 
prepare, evolve, and read out the quantum state of the circuit^i^i^i^. One of the most successful readout methods 
Q ■ used this far is based on the switching of a Josephson junction from the zero voltage state to the finite voltage state. 
O In this letter we analyze this switching process in some detail, focusing on the reliability of the detection method in 
. the presence of external noise. We show that the existence of an unstable limiting cycle in the phase space dynamics 
' of the Josephson junction leads to late retrapping (long measurement time) and false switching events, and that such 
^ , an unstable cycle can be avoided with proper choice of parameters, while maintaining the desired overdamped phase 
• dynamics. 

The switching current detector works on the principle that the quantum states one wants to differentiate have 
different critical currents lon — max[dEn/d(j)], where E„ is the energy of eigenstate n, and (p is the external phase 
of the circuit, which can be changed by application of an external current. Restricting ourselves to the two lowest 
energy states of the circuit (a qubit) we quickly ramp the current to a value Ip, where /qo > Ip > Iqi- The circuit 
' will evolve differently, depending on whether it is in state or in state 1. If the qubit is in state 0, the circuit will 
, not switch, meaning that it remains in the zero voltage state where the phase is trapped in a local minimum of the 
' potential E'o (</>)■ If the qubit is in state 1, it will switch, meaning that the phase escapes from this minimum, and 
^ [ evolves to a "free running" state with finite voltage. Reading the absence or presence of this voltage constitutes a 
I ' determination of the quantum state. 
' O I A good switching detector should have a bistability, characterized by two stable attractors in the phase space of the 
^ ' non-linear circuit. As we shall explain, such a bistability makes possible a "sample and hold" approach to quantum 
measurement^. For the Josephson junction qubit detectors studied thus far, bistability has been realized by switching 
Josephson junctions with underdamped phase dynamics at all frequencies. A junction with overdamped dynamics 
. !^ ' frequencies would not have a bistability. Overdamped dynamics is however desirable from the point of view 

, of minimal phase fluctuations (which cause dephasing of the qubit) and rapid reset time of the switching detector. 
5_j Therefore, an interesting model in the context of qubit detectors is the Josephson junction shunted by a series RC 
' circuit, which has the property of overdamped phase dynamics at high frequencies, and underdamped dynamics at low 
frequencies. This model has been studied by several authorsSiSiifliiLiS. Kautz and Martinis showed that bistability is 
possible in this model, with one state characterized by "phase diffusion" where the phase is thermally kicked between 
minima of the potential, and the other state is that of a free running phase. 

A good detector has high readout fidelity and short measurement time. To achieve high fidelity the circuit should 
be insensitive to all sources of noise which may cause false switching events, so that the switching probability is only 
determined by the quantum probabilities. Short measurement time is necessary so that the measurement occurs before 
the qubit has time to change state (relaxation or excitation). In the switching detectors studied thus far, the fidelity 
and speed is thought to be limited by macroscopic quantum tunneling of the phase in the Josephson potentiaU^*i^. 
Here we consider switching of junctions which are very overdamped at high frequencies, with small capacitance and 
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FIG. 1; The model circuit (a) and the apphed current pulse (b). 



shunting impedance much less than the quantum resistance {h/e^). In this case a classical description of the non- linear 
phase dynamics is adequate. 

The most simple qubit circuit is a Cooper pair transistor. However, in order to study the switching process it is 
sufficient to consider only a single small capacitance Josephson junction. The Josephson junction has critical current 
/o, capacitance Co and is biased by a current source I{t). In parallel with the junction is a resistor i?i and a series 
combination of a resistor R2 and a capacitor C2. See Fig.nja). 

The classical dynamics of our system can be described in terms of cj) (the phase over the junction), v and vq (the 
ratio of the voltages over the junction and over C2 to RiIq). However, for the small junctions considered here, the 
current through Cq can to a very good approximation be neglected and we are left with 
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are Gaussian distributed and obey (ini,2(i'i)*ni,2(i2)) = 
= i?i ■\/2e7oC'o7^ is the low frequency quality factor, Qi — 
R2^)~^ \/2eIoCo/h is the high frequency quality factor and p — RiCq/R2C2- The combination {i 
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the force associated with the well known tilted washboard potential 

Fig.|21shows the {(j3,vc) phase portrait for a finite bias current ihoid = 0.80 which is applied during readout. Two 
cases are shown, both having Qo = 3.5, with (a) Qi = 0.02 and p = 0.014 and (b) Qi = 0.04 and p = 0.006. In both 
cases there are stable fix points A (located at (arcsin(i^o;(i) + 2im, 0) corresponding to local minima of the washboard 
potential; shown with filled red circles), and saddle points S (located at (tt — arcsin(i;io/(i) + 27m, 0) corresponding 
to local maxima of the washboard potential; shown with red/black cirlces). There is also a stable limiting cycle 
B (located close to vc ~ ihou, not shown in the figures). The regions of dotted red trajectories are the basins of 
attraction for the different A's (non-switching events), the region of solid blue trajectories is the basin of attraction 
for B (switching events). The boundary between the switching and the non-switching regions are shown with thick 
black trajectories. The arrows point in the direction of positive time evolution. Figs. EJa) and (b) show a striking 
difference in topology: in (b) the boundary is an unstable limiting cycle C (with Q < vc < Vsw being periodic in 0), 
whereas in (a) there are several boundaries (all ending in the different S's) which separate the different red regions. As 
Ql is decreased, C approaches S. When C collides with S and vanishes, the black trajectory flips origin to vc = —00, 
and the blue basin stretches down and separates the different red basins as shown in (a). For the parameters used 
here, this flip occurs close to Qic ~ 0.029. 

The sample and hold switching measurement is realized by a fast initial pulse /p, with duration as shown in Fig. 
1(b). After this pulse, the phase space coordinate {(l){Tp),vc{Tp)) will lie close to the boundary between the basins of 
attraction. The existence of the unstable boundary C imples two things: flrst, a long time is required for trajectories 
near the boundary to separate into the two basins of attraction. Second, events can jump between consecutive (red) 
non-switching basins without passing the blue switching region. When Qi increases, both the separation time and the 
distance between A and C will increase. Hence for Qi > Qic (Fig.^fb)), the topology gives rise to "late retrappings" 
and "late switches" which increases the measurement time and result in errors. On the other hand, for Qi < Qic 
(Fig. 121a)) the boundary is a saddle point flow, and there is no C, and consequently no late retrapping or late switching. 
At finite temperatures late events are still possible, although the probability of occurrence is much smaller than for 
Ql > Qic and decreases with smaller Qi. 
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FIG. 2: (Color online) (0,wc) phase space flows for Qo = 3.5, hou = 0.80 and (a) Qi = 0.02, p = 0.014; (b) Qi = 0.04, 
p = 0.006. 
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FIG. 3: (Color online) Simulated {(j>,vc) values at t = 6.5 ns for 1000 events, ip giving 50% switching, Qo — 3.5, ihoid ~ 0.80 
and Tn = 30 mK. Switching (non-switching) events are marked with blue plus-signs (red dots), miscounts with diamonds. In 
(a) Qi = 0.02, p = 0.014 and in (b) Qi = 0.04, p = 0.006. 



To demonstrate this sensitivity to phase space topology we have carried out numerical simulations with experi- 
mentally relevant parameters; Tp = 5 ns, Iq — 100 nA, Cq — 10 fF, T„ — 30 mK, and a 2% random variation in 
the value of Ip which gives 50% switching. Initially, ((/)(0), wc(0)) = (0,0) which corresponds to the phase being 
trapped in a minimum of the washboard. Adding the current pulse and hold and thermal noise will evolve the system 
to a distribution {{(j)(t),vcit))}. The speed of the detector is determined by the time at which the distribution 
distinctly separates into the two basins of attraction. In the optimal case, = Tp. 

In Fig.|3Ia,b) we have plotted (</',i'c) distributions for 1000 events at i = 6.5 ns. Switching (non-switching) events 
are marked with blue plus-signs (red dots), events that change basin after 6.5 ns (i.e. miscounts) are marked with 
diamonds in the respective color. We see that the distribution for Qi = 0.02 (topology of Fig.|2fa)) is clearly separated 
into the two basins of attraction, with only a handfuU of events close to the boundary, and no miscounts. Note that 
all non-switching events arc in the first minimum, and all switching events emanate from the first saddle point, with 
no late switches or late retrappings. On the other hand, for Qi = 0.04, many events lie close to C and there are 
several late switches and late retrappings. There are 34 miscounts, and for higher Qi there will be many more such 
events. 
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In summary, we have shown that the phase space topology plays a crucial rule in the fidelity and speed of switching 
current readout for quantum bits. We find that it is possible to achieve almost optimal speed and very high fidelity 
with overdamped phase dynamics at high frequencies. 

We would like to acknowledge helpful conversations with D. Vion and S. Ostlund. This work was partially supported 
by the EU project SQUBIT-2. 
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